Constant Flux Relation for diffusion limited cluster— cluster aggregation 



OO 
O 

o 

(N 
C 

o 

S 1 

o 

£ 

a 

a 

c 

o 
o 



> 

m 
cn 

o 

oo 
o 



X 



Colm Connaughton, 1,2 ^ R. Rajesh, 3 '0 and Oleg Zaboronski 2, [l 

1 Centre for Complexity Science, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK 
2 Mathematics Institute, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK 
3 Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai-600 113, India 

(Dated: June 20, 2008) 

In a non-equilibrium system, a Constant Flux Relation (CFR) expresses the fact that a constant 
flux of a conserved quantity exactly determines the scaling of the particular correlation function 
linked to the flux of that conserved quantity. This is true regardless of whether mean-field theory 
is applicable or not. We focus on cluster-cluster aggregation and discuss the consequences of mass 
conservation for the steady state of aggregation models with a monomer source in the diffusion- 
limited regime. We derive the CFR for the flux-carrying correlation function for binary aggregation 
with a general scale-invariant kernel and show that this exponent is unique. It is independent of 
both the dimension and of the details of the spatial transport mechanism, a property which is very 
atypical in the diffusion-limited regime. We then discuss in detail the "locality criterion" which must 
be satisfied in order for the CFR scaling to be realisable. Locality may be checked explicitly for the 
mean-field Smoluchowski equation. We show that if it is satisfied at the mean-field level, it remains 
true over some finite range as one perturbatively decreases the dimension of the system below 
the critical dimension, d c — 2, entering the fluctuation-dominated regime. We turn to numerical 
simulations to verify locality for a range of systems in one dimension which are, presumably, beyond 
the perturbative regime. Finally, we illustrate how the CFR scaling may break down as a result of 
a violation of locality or as a result of finite size effects and discuss the extent to which the results 
apply to higher order aggregation processes. 

PACS numbers: 05.20.-y, 47.57.eb, 68.43.Jk, 61.43.Hv 



I. INTRODUCTION AND MOTIVATION 

Consider a collection of particles undergoing some spa- 
tial transport process which, upon encountering each 
other, coalesce irreversiby with some probability. Such 
a situation arises in a great variety of seemingly unre- 
lated branches of science (see [l| for an overview). Some 
of the most obvious examples are found in astrophysics 
0, aerosol physics [H and polymer chemistry [j]. Less 
obvious examples arise from granular media [f| , the struc- 
ture of drainage networks [6|, |7| and sandpile models of 
self-organised criticality This diverse range of appli- 
cations is one reason why models of systems of diffusing 
particles which aggregate upon contact have been exten- 
sively studied since the seminal work of Smoluchowski 
laid the foundations for their analysis. A second reason 
for the enduring interest shown by the scientific commu- 
nity in aggregation models is that they provide simple 
examples of a surprising array of non-trivial phenomena 
in non-equilibrium statistical mechanics making them an 
attractive theoretical proving ground. 

Two situations are commonly encountered, depending 
on the application. One may start with a specified initial 
distribution of cluster sizes and study how it decays in 
time. This is sometimes referred to as free aggregation. 
Alternatively one may start with an empty system and 
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add monomers at a given rate. This is called aggregation 
with a source. Due to irreversibility of the coagulation 
process, free aggregation is an entirely dynamic problem 
with no stationary state. On the other hand, aggrega- 
tion with a source may produce a stationary distribution 
of particle sizes in the limit of large time. Stationarity 
comes about as follows: the rate of decrease of the den- 
sity of clusters of a given size via coagulation to form 
larger ones is balanced by the generation of clusters of 
that size via coagulation of smaller ones. Such a balance 
is possible only because the source continually replen- 
ishes the available pool of small clusters. Clearly such 
a stationary state is not an equilibrium state since there 
is no detailed balance. Rather it is a flux state charac- 
terised by a constant flux of mass through the space of 
cluster sizes. On a technical note, since both diffusion 
and aggregation conserve total mass, the constant influx 
of monomers results in a linear increase in the average 
mass. While this driving occurs at the smallest mass in 
the problem, the aggregation process transfers this mass 
to larger and larger mass scales. Thus, strictly speaking, 
such systems are quasi-stationary at large times: small 
masses reach a stationary distribution but time-evolution 
proceeds indefinitely at the largest masses. To attain a 
truly stationary state, one should introduce a cut-off at 
some large cluster size above which clusters are removed 
from the system. In this paper, we concern ourselves 
exclusively with aggregation problems with a source. 

The most basic quantity of interest is the average mass 
density, (N(m,x,t)), which tells us the average number 
of clusters of a given mass, m. Typically, a system of ag- 
gregating particles exhibits two regimes of behaviour as 
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a function of the spatial dimension. A critical dimension, 
normally two for systems undergoing diffusive transport, 
separates these regimes. In higher dimensions, the dy- 
namics is typically reaction limited and a mean-field de- 
scription is appropriate. This mean-field description is 
given by the Smoluchowski kinetic equation which de- 
scribes the time evolution of N(m, x, i). In lower dimen- 
sions, the dynamics is typically diffusion limited. Diffu- 
sive fluctuations are strong and a mean-field description 
is no longer possible. A huge amount is known about 
the average mass density in the mean-field case [|| from 
exact analyses [lfj| and extensive numerical simulations 
of the Smoluchowski equation. Relatively less is known 
about the mass density in the diffusion limited regime but 
several models have been solved exactl y o r treated ap- 
proximately by field-theoretic methods [ill . [l2| . Almost 
nothing is known about higher order correlation functions 
in the diffusion-limited regime, despite the fact that they 
encode the details of the fluctuations which dominate the 
dynamics. This paper concerns itself with such higher or- 
der correlation functions, albeit some rather special ones. 

The special correlation functions which we consider, 
can be referred to as the flux-carrying correlation func- 
tions. For a given aggregation model with source which 
attains a constant-flux stationary state at large times, 
there is a particular correlation function associated with 
mass transfer. In the turbulence literature, where trans- 
fer of energy is analogous to transfer of mass, it is very 
well known that constancy of the energy flux determines 
exactly the scaling of the flux-carrying correlation func- 
tion (see chap. 6 of [l3[ for example). This fact is the 
basis for the Kolmogorov 4/5-Law for three-dimensional 
hydrodynamic turbulence and its various incarnations in 
other turbulent systems. While the 4/5-Law has become 
central to the modern understanding of turbulence, the 
fact that a similar exact result is available for other non- 
equilibrium systems, in particular for aggregation sys- 
tems, has hardly been taken advantage of. The pur- 
pose of the present article is to address this issue. In 
previous work [l4| . we showed how a conservation law 
leads to an exact scaling exponent for the flux-carrying 
correlation function for a broad class of non-equilibrium 
systems which included aggregation, referring to such a 
constraint as a "Constant Flux Relation" (CFR). In the 
present article we focus entirely on the consequences of 
CFR for aggregating particle systems, leaving the origi- 
nal hydrodynamic analogy behind. In the process of ver- 
ifying the CFR for a broad set of aggregation models in 
the diffusion-limited regimes we will present a number 
of somewhat counter-intuitive numerical results which 
would be very difficult to understand without any prior 
understanding of the CFR. 

The layout of the paper is as follows. We first define 
the model and give a heuristic derivation of the CFR 
scaling (Sec.|TI|. We then provide an accurate derivation 
fSec. HIl| ) which makes explicit the assumptions involved, 
in particular the assumption of locality which we then 
discuss in detail (Sec. HV)) . Sec. |V]then reports the re- 



sults of a large number of numerical simulations which 
verify the CFR scaling for a range of aggregation ker- 
nels, expose finite size effects, test the locality condition 
in the diffusion limited regime in one dimension where an 
analytic approach is lacking and demonstrate the lack of 
dependence of the CFR scaling on the details of the dif- 
fusion. Finally we extend the discussion to higher order 
aggregation processes (Sec. I VI I) . We close with a brief 
summary of the results. 



II. MODEL DEFINITION AND HEURISTIC 
CFR 

Consider a d-dimensional hypercubic lattice occupied 
by point size particles carrying a positive mass. Multi- 
ple occupancy of a site is allowed. Given a certain con- 
figuration, the system evolves in time via the following 
processes. 

• Diffusion: A particle hops with a mass dependent 
diffusion rate D(m) to a randomly chosen nearest 
neighbour. 

• Coagulation: Two particles of masses mi and m.2 
on the same lattice site coagulate at rate A (mi, TO2) 
to form a particle of mass mi + m-x- 

• Input: Particles of mass mo are injected at rate 
J/mo uniformly and independently in space. 

The initial condition is one where the lattice is empty. 
We shall call this model the mass model (MM). 

We will restrict ourselves to the case where the re- 
action rate A(mi,m2) is a homogeneous function of its 
arguments, i.e., 

A(Ami, Am 2 ) = A /3 A(mi,m 2 ), (1) 

where f3 is the homogeneity exponent. The diffusion con- 
stant, D(m), will be assumed to have the property 

D ( m ) = ( m_Y (2] 
D{m ) \m J { ' 

Thus, in addition to the different rates, the model has 2 
parameters: the homogeneity exponent /3 and the diffu- 
sion exponent k. In the large time limit, as described in 
the introduction, this model tends to a statistically sta- 
tionary state characterised by a constant average flux of 
mass from small clusters to large ones. 

In [l4| we presented quite a general argument to deter- 
mine the scaling of the flux-carrying correlation function 
for a broad class of non-equilibrium systems which reach 
a constant flux stationary state. In this article, in the 
interest of clarity, we will briefly review the argument 
heuristically for the specific case of particle aggregation. 

Schematically (we shall write down an accurate expres- 
sion in Sec. IIIip . the transfer of mass between coalescing 
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clusters is described by an equation of the form: 

9 i at u\\ dJm 
di {mNmit)) = W 



dmidm,2 m A(mi, m 2 ) C(mi, m 2 ) c>o;i,2, (3) 



where 5o-i,2 is shorthand notation for S(m — mi — m 2 ). 
The right hand side defines the mass flux, J m , in the 
space of cluster sizes. C(mi,m 2 ) is proportional to the 
probability of having two clusters with masses mi and m 2 
meet at the same point in space. This is the flux-carrying 
correlation function since it mediates the transfer of mass 
in the system. Note that the flux-carrying correlation 
function is not an esoteric object. It has a clear and 
intuitive physical meaning 

In the statistically stationary state, dN g t ^ = so that 
J rn is a constant, independent of m. Simply counting 
powers of m would then lead us to expect that 



C(mi, 7712) ~ m 



-0-3 



(4) 



This heuristic scaling argument is the CFR at the most 
basic level: mass conservation fixes the scaling of the flux- 
carrying correlation. The remainder of the paper will be 
devoted to making this heuristic argument precise and 
identifying its limitations. 



III. IMPROVING ON THE HEURISTIC CFR 

In this section we arrive at Eq.(|3]) more carefully. 
Starting from the lattice model, it is relatively straight- 
forward to write down the evolution equation for the 
different correlation functions. A full exposition can be 
found in [l2j]. Skipping the details, we write directly the 
equation for (N(m, x, t)}, the average number of particles 
of mass m at position x at time t: 



d \ I 

— -D(m)V 2 ){N(m)) = — 5(m-m ) (5) 

at J mo 



dmidm 2 X(mi,m 2 ) C(mi, m 2 ) <5o ; i,2 

1 

dm\dmi A(m, mi) C(m, mi) ^2;0l 



dmidm 2 A(m 2 , m) C(m 2 , m) ^i ;2 o- 



For simplicity, we suppress x and t dependences and 
adopt the reduced notation for the 5— functions defined 
after Eq.Q. C(mi,m 2 ), the flux-carrying correlation 
function, is defined: 

C(mi,m 2 ) = (N(m 1 ,x 1 t)N(m 2 ,x,t)) (6) 
~~ A-z^ 7 ™ 1 ~~ m i)(N(mi,x, t)), 

Ax being the lattice spacing. Let us explain the terms 
in Eq (O one by one. 



The V 2 term accounts for particle diffusion which may 
be mass dependent. For spatially homogeneous statistics, 
this term is zero. The first term on the right hand side ac- 
counts for influx of particles of mass ttiq. The remaining 
terms account for aggregation processes. To explain the 
meaning of C(mi,m 2 ), we first consider how it relates 
to the mean-field Smoluchowski equation. Mean-field 
theory requires two assumptions. Firstly, correlations 
are absent so we may write {N(mi,x, t)N(m^, x, t)) as a 
simple product of densities, (N (mi)) (N (1712)) ■ Secondly, 
densities are high so we may neglect the (7V(mi)) term 
relative to (N(mi))(N(m 2 )). In the diffusion-limited 
regime, C(mi,m 2 ) has an important probabilistic inter- 
pretation. Writing the averaging process explicitly: 

00 

C(mi,m 2 ) = ^2 F(N(m 1 ,x) = Nx,N(m 2 ,x*) = N 2 ) 

Nx,N 2 =l 

X {NiN 2 ~ S mum2 Ni) . 

This is the average number of pairs of particles with 
masses mi and m 2 on a site, with the delta function ac- 
counting for double counting of particles of equal mass. 
In the low density (diffusion-limited) regime, 

C(mi,m 2 ) w ¥(N(rnx,x) = 1, N(m 2 , x) = 1), (7) 

the probability that two particles of masses mi and 
m 2 meet at a site. Thus the flux-carrying correlation 
function is not an esoteric object and has a very nat- 
ural physical meaning. Having understood the mean- 
ing of C(mi,m 2 ), the second term on the right-hand 
side of Eq. ([5]) accounts for the creation of particles of 
mass m at x through aggregation of 2 particles at x. 
The third and fourth terms account for the decrease 
of N(m,x,t) through aggregation with other particles. 
These latter two terms are identical under relabeling 
(mi,m 2 ) —> (m 2 ,mi) and are usually written as a sin- 
gle term. We write them this way for reasons which will 
become obvious below. 

To simplify the equations, we introduce 7(mi,m 2 ;m) 
defined as: 

I(mi,m 2 ;m) = A(mi, m 2 ) C(mx, m 2 ) 5 0; i,2 (8) 

As already mentioned, in Eq. ([5]) the diffusion term 
drops out by spatial homogeneity. Then, for m > mo we 
can write Eq. ([5]) as 



d(N(m)) 



dmidm 2 



dt jo 
-I(m 2 , m; mi) — I(m, mi; m 2 



I (mi, m 2 ; m) 



(9) 



In the steady state, we set the left hand side to zero. 
To solve this equation, we need to balance out the plus 
and the minus terms on the right hand. As written, it 
is difficult to see what the solution is because each term 
comes with a different delta function. The balance can 
be made explicit by the Zakaharov transform (ZT) [l5l |. 
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Leave the first term as it is. Make the following trans- 
formation of the second integral 



1TL2 



mmi 

m 2 

o 

m 
m 2 



(10) 
(11) 



The Jacobian of the transformation is (m/m 2 ) 3 . Perform 
the analogous transformation of the third integral (see 
[IB]). Now look for homogeneous solutions, i.e., 

C(Ami,Am 2 ) = K h C(m 1 ,m 2 ) 



(12) 



Using this and the homogeneity of A, we obtain 



= / dm\dm, 2 I{m\,m,2\Tri) (m y — m\ — m^) (13) 
Jo 

where y = —h — (3 — 2. 

Due to the delta function in Eq©, I is non-zero only 
when mi + m 2 — m. If the term in the square bracket 
is zero when / is non-zero, then the equation is satisfied. 
Thus, y = 1 is a solution. This implies 



-/3-3. 



(14) 



It can be easily shown that this is the unique homeoge- 
neous stationary solution of Eq. |9]). Introducing rescaled 
variables, x\ — mi/m and x 2 — m 2 /m and using the as- 
sumed homgeneity of C(mi, 7712), Eq. (fTB")) can be rewrit- 
ten as 

= m 1+h+l3+y f dx 1 dx 2 I{x ll x 2 ]l){l-x v l -x v 2 ) (15) 



Due to the delta function in I(x\,x 2 \ 1), the integrand 
is zero unless x\ + x 2 = 1. When x\ + x 2 = 1, the 
integrand clearly vanishes for y = 1. To show that 
this is the only value of y for which the integral is 
zero, we show that for y ^ the integrand is sign def- 
inite on the domain of integration so that the integral 
is not zero. From the definition, Eq. ([8]), I(xi,x 2 ;l) 
is clearly positive. It remains to consider the function 
f(xi,x 2 ) = \ — x\—x\. For y > 1 the fact that Xi E (0, 1) 
implies that x\ < 1 so that x\ +x 2 < x\ + x 2 = 1. Thus 
f(xi,x 2 ) > and the integrand is everywhere positive. 
Likewise, for y < 1, € (0,1) implies that x\ > 1 so 
that x\ + x\ > x\ + x 2 = 1. Thus f(xi, x 2 ) < and the 
integrand is everywhere negative. Thus, for y ^ 1 the 
integral does not vanish and the only solution is y = 1. 

One may make a curious observation: the diffusion 
constant does not play any role. This is counter to the 
usual intuition in reaction diffusion systems which holds 
that diffusion is unimportant for dimensions greater than 
upper critical dimension and all-important for dimen- 
sions lower. Here, we have shown that the 2-point cor- 
relation function is independent of dimension and of the 
spatial transport mechanism. 

It must be pointed out that these manipulations are 
correct provided each of the integrals in the evolution 
equation are convergent. This condition referred as the 
locality condition has to checked separately. This will be 
discussed next. 




FIG. 1: The support of the integrand of Eq. ([5]) in the mi, 7712) 
plane. 



IV. LOCALITY: WHEN IS CFR REALISABLE? 

To obtain the formal scaling solution [Eq. (fT4")) ] for 
the evolution equation (Eq. (O), some implicit assump- 
tions were made. These assumptions will be referred to 
as locality condition, the terminology being borrowed 
from wave turbulence. Unless, these assumptions can 
be proved or checked numerically, the scaling solution 
should not be expected to hold. In this section, we ex- 
plain the locality condition in detail. 

For the scaling solution with exponent given by 
Eq. (fT4"|) to be physically realisable, it must yield a con- 
vergent integrand on the right hand side of Eq. ([5]) , before 
any changes of integration order are made. Otherwise, 
divergences cancel leaving a finite contribution. 

To study this, let us write the two point function as 



C(toi,to 2 ) = (mim 2 ) h/2 (t> 



rrii 
m 2 



(16) 



thus introducing the dimcnsionlcss scaling function 4>{x) . 
(f>(x) has the symmetry property 4>{x) = 0(x _1 ). To 
check convergence, it is not enough to know just the 
degrees of homogeneity but rather we require to know 
limiting behaviour of various quantities in the integrand. 
Suppose 



A(mi,m 2 ) 
4>{x) 



i 2 , for m 2 » mi, (17) 
for x < 1. (18) 



The exponents /j and v are determined by the model 
under consideration and must satisfy [i + v — [3. The 
behaviour of the scaling function (j){x) as x — > 0, as de- 
termined by the exponent a, is something which we do 
not apriori know. 

The support of the integrand in Eq. ([5]) is shown in 
Fig. [TJ We may integrate once and consider the integral 
as an integral in mi only. By scale invariance, we need 
to check convergence only at the endpoints of the range 
of integration. The analysis was done in Ref. [l6[ in the 
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mean field limit. Following the analysis of [lq , as mi 
oo, the behaviour of the integrand is given by 



A(to, mi)(mmi) h ^ 2 cf> 



s mi / \mi 

The integral is convergent at infinity if 

-h/2>v+l-a. (19) 

For the behaviour at m\ — > 0, there is a cancellation of 
leading order terms : 

A(mi, to — toi)C(toi, to — toi) — A(toi, to)C(toi, to) 

~ i x ( m i, x ) c ( m i, x )] x = m + o(m\), 



TOi 



(9.T 



^(mia^/W^Vl +o{m\). 

\ X ' -I a;— m 



The integral is convergent at if 

- h/2 < 2 + fi + a. 



(20) 



Putting together Eq. ([19]) and Eq. ([20]) . a convergent col- 
lision integral requires that the interval [v+l—a, fi+2+a] 
should have positive width. The width of this interval is 
2<T + /i — v + 1 . Thus a convergence requires 



!)• 



(21) 



It is easy to show that if this interval exists, the expo- 
nent —h/2 lies within it assuring the validity of the CFR 
solution. 

At the level of mean field theory a = since 
C(toi,to,2) is simply proportional to the product of the 
1-point densities. This case was worked out in detail in 
Ref. [lj| and is consistent with Eq. (|2"Tj) . 

Thus the rigorous verification of CFR in MM requires 
the knowledge of the small- X behaviour of the scaling 
function. The latter can be often studied using pertur- 
bative methods. For instance, consider constant kernel 
MM, [i = v = 0. If dimension of the physical space 
is two, mean field approximation is applicable, perhaps 
modulo logarithmic corrections. Within mean field ap- 
proximation, (7 = and criterion Eq. (|21[) is satisfied. 
Hence, CFR holds in two dimensions, meaning in this 
case, C(toi,to,2) ~ to~ 3 , and logarithmic corrections are 
absent. Consider now constant kernel MM in d = 2 — e, 
where e > 0. The dynamics of the model is governed 
now by a fixed point of renormalization group. The or- 
der of the fixed point is e. Scaling exponents can be 
now computed using e-expansion. As er(e = 0) = 0, 
cr(e) = die + 0(e 2 ), where o\ is a constant. Assuming 
that (7i 7^ 0, one can re-write locality criterion as 



eai > -\, 



(22) 



which is satisfied, if e is small enough. In this case all 
e-corrections to h = —3 vanish and CFR holds. 




FIG. 2: The variation of n2(m) with m is for = —0.5 (top 
curve), —0.25,-0.10,0 (bottom curve). The solid lines cor- 
respond to the exponent as predicted by CFR [see Eq. (|14p ]. 
Curves have been slightly shifted for clarity. 



This perturbative argument demonstrates that sys- 
tems which are local at the mean-field level remain so for 
some time as one decreases the physical dimension into 
the diffusion-limited regime. It does not tell us much 
about whether locality holds by the time one reaches the 
next physically relevant, integer valued dimension below 
the critical dimension since this is presumably beyond 
the perturbative regime. In general, this is a difficult 
problem. The only case which we are aware of which can 
be handled analytically is the constant kernel ((3 = 0). 
It may be shown to be true for d > 2 where mean field 
holds [l(| , in d = 2 due to a cancellation of logarithmic 
corrections [l2| and in d < 2 by an exact solution [l8[ . 

Lacking an analytic approach for other kernels, one 
must rely on numerical simulations to measure the value 
of a for particular systems. We perform a systematic 
numerical investigation of locality in one dimension for 
several kernels in Sec. IV CI 



V. NUMERICAL SIMULATIONS OF CFR 

A. Numerical measurements of CFR exponent 

In this section, we present results of numerical simula- 
tions directly measuring the exponent h given in Eq. (|14[) 
It is in one dimension that the effects of fluctuations are 
the strongest. Hence, if the mean field scaling predicted 
by CFR is violated, then it will be violated in one dimen- 
sion too. For this reason, all the numerical results that we 
show will be monte carlo simulations for one dimensional 
lattices. 

In our simulations we investigated the following repre- 
sentative kernels: 



A(toi, TO2) = rrii + 



TO 



(23) 
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FIG. 3: The variation of n2(m) with m for /3 — 0.50 (top 
curve), 0.75, 1.00, 1.50 (bottom curve). The solid lines corre- 
spond to the exponent as predicted by CFR [see Eq. (|14p ]. 
The dotted line corresponds to exponent for /3 = 0. Curves 
have been slightly shifted for clarity. 



FIG. 4: The lattice size dependence of H2{rn,L) is shown 
for n = 3 and /3 = 1.5. The different lattice sizes are 
L = 100, 1000, 10000. The crossover moves to the right with 
increasing L. Inset: The curves are scaled as in Eq. (|27p . The 
straight lines correspond to slopes —7/2 and —2. 



A(mi,m 2 ) = {m im2 fl 2 , (24) 
A(mi,m2) = max[mi, m2] ly min[mi, 7712]^, (25) 

which we shall refer to as the additive, multiplicative and 
mixed kernels respectively. As far as the value of h was 
concerned, the results were identical for all three kernels. 
Hence, unless stated otherwise, our figures present results 
only for one of them, namely the additive kernel, Eq (|2U)) . 

What is convenient to measure in simulations is not 
(N(nii)N(m2)), but the quantity, 

/>OG 

7r 2 (m) = / dm 1 {N{m)N(m 1 )). (26) 

J m 

CFR predicts that 7r2(m) scales as 7T2(m) ~ m~ 2 ~^ . 

In Fig. [51 the variation of 7T2 (to) with to is shown for 
/3 = 0, -0.10, -0, 25, -0, 50. The solid lines are the CFR 
results. As can be seen, there is excellent agreement, 
confirming that CFR holds when j3 < 0. Fig. [31 shows the 
variation of 7t 2 (to) with to for (3 = 0.50,0.75,1.00,1.50. 
The solid lines are the CFR results. The CFR exponent 
is obtained for small and intermediate masses but there 
is a clear cross-over to another behaviour at large masses. 
This can be understood as a finite size effect which we 
shall discuss in section IVBl 

At this point it is appropriate to make some comments 
on the correspondence with mean field theory. For the 
additive kernel, Eq ([23j) . at the mean field level it is be- 
lieved [l!| that = 1.0 corresponds to the threshold for 
instantaneous gelation so that analytical understanding 
of the solutions of the Smoluchowski equation for (3 > 1 is 
very difficult. Notwithstanding the cross-over to another 
regime at large masses, it is very interesting that the 
CFR exponent is observed over some considerable range 
even for j3 > 1. To the best of our knowledge, nothing 



is known about the behaviour of gelling kernels in the 
diffusion limited regime or in the presence of a monomer 
source. In this light, the results of Fig. [3J pose many in- 
teresting questions such as whether there is any remnant 
in the diffusion-limited regime of the catastrophic singu- 
larity which occurs in the mean-field equation at /3 = 1. 



B. Finite size effects 

Let us consider why the lattice size should affect the 
CFR scaling. At this point it is useful to recall two things. 
Firstly, the recurrence property of random walks plays a 
crucial role in determining the statistics of aggregation 
in the diffusion limited regime. Due to recurrence, heavy 
particles develop "zones of exclusion" around them as 
they grow, resulting in strong anti-correlations between 
heavy particles. Secondly, recall that the CFR exponent 
quantifies the decreasing probability of two heavy parti- 
cles meeting each other. Due to the presence of zones 
of exclusion, this probability decreases faster for heavy 
particles than the product of one-point densities would 
suggest. 

Although zones of exclusion grow larger as particles 
get heavier, in an infinite system there are always enough 
heavy particles to maintain the CFR scaling over all mass 
scales. In a finite system, however, these zones of exclu- 
sion become limited by the system size eventually. Once 
this happens, heavy particles start to meet each other 
more often than would be expected from CFR since thay 
can no longer grow their zones of exclusion any larger. 
Thus a finite size cross-over occurs and results in a shal- 
lower scaling as evident from Fig. [3J No such crossover 
occurs for the j3 < cases shown in Fig. [5] since for 
f3 < large mass clusters become progressively less reac- 
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\ v 

p\ 


-0.250 


-0.125 


0.000 


0.125 


0.250 


0.375 


0.500 


-0.25 


1.33 


1.33 


1.35 


1.44 


1.57 


1.69 


1.82 


0.00 


1.34 


1.32 


1.34 


1.45 


1.58 


1.70 


1.83 


0.25 


1.32 


1.31 


1.34 


1.46 


1.58 


1.70 


1.83 


0.50 


1.31 


1.31 


1.33 


1.46 


1.59 


1.70 


1.82 



TABLE I: The numerical values of a — h/2 are shown 
for different v and (3. The kernel used in A(mi,m2) = 
max(mi, 7712)" min(mi, m2) M . The errors in the values are 
±0.02. 



tive which acts to counterbalance the growth of zones of 
exclusion due to recurrence. 

The argument above does not explain why finite size ef- 
fects should lead to a the scaling corresponding to (3 = 
indicated in Fig. [3] We suggest the following heuristic 
argument. For a finite system size, C{mi,m,2) for 'large' 
masses is contributed to by configurations consisting of 
two heavy particles which have been in the system for 
times ^> L 2 , so that they are strongly-anticorrelated. 
Hence these two particles effectively interact with each 
other at infinite rate, with effective diffusive jumps of 
the size equal to system size. Hence C(mi,TO2) behaves 
as if beta=0 at these masses. Since the mass flux is car- 
ried by the meetings of these super-heavy particles, it is 
presumably highly intermittent. It is then intuitive that 
the constant flux argument should fail to describe this 
regime. 

Given that we expect to see CFR scaling for small 
masses and (3 — behaviour for large masses, we expect 
that 7T2(to, L) should have the form: 



7r 2 (m, L) 



1 



L1+2//3 



m 



(27) 



-2-13 



where the scaling function f{x) varies as f(x) < 
when x — > and f(x) ~ m~ 2 when x 3> 1. The crossover 
mass m c is given by ~ L, or m c ~ L 1 ^ . 

In Fig. [4] we study the variation of ^(m, L) with m 
for fixed (3 and different L. The (3 value is chosen to be 
/3 = 1.5. As expected from the preceeding discussion, 
the crossover point moves to the right with increasing 
L. In the inset, the data is scaled according to Eq. ([27]) 
and excellent collapse is obtained. The large and small x 
behaviour of the scaling function behaves as predicted. 



C. Numerical validation of locality criterion 



As stressed in Sec. IIV1 aside from a couple of special 
cases we do not know whether the locality criterion is 
satisfied in one dimension or not. We now present nu- 
merical measurements of the exponent a in Eq. (|17p to 
address this issue. We choose our lattice size sufficiently 
large to avoid any question of the finite size effects dis- 
cussed in the previous section influencing the exponents. 



a 



10* 
10° 
10" 2 

10- 4 

10" 6 
10" 8 



10 



10 



10 



12 




10" 



P = 0; K= 0.0 
p = 0;k = -0.5 

10 1 



lo- 



ur 1 



FIG. 5: The dependence on diffusion constant of of iT2(m,L) 
is shown for n = 3 and f3 = 0. The bottom curve corresponds 
to D(m) ~ m° [k = 0] and the top curve corresponds to 
D(m) ~ m -1,/2 [k = — 1/2]. The curves have same slope. 
Inset: The variation of {N(m)) with m is shown. There is a 
strong dependence on ft. The solid lines have slope (4 + ft)/3. 



We use the mixed kernel, Eq. ([25]) so as to be able to vary 
v and /1 independently. What is measured numerically 
is (N(mi)N(m2)) when mi is kept fixed and rri2 ^> vn\. 
Then (N(mi)N(m2)) ~ m^ 2+c 'm^ 2 ° . In our simula- 
tions we keep m\ fixed at 5mo and take large and 
measure a — h/2. The results of a systematic set of nu- 
merical experiments are shown in Table |TJ What one sees 
is that a — h/2 is independent of (3 and dependent only 
on v. The numerics suggest that: 



^4 r . 

cr - - = - +max[;y,(Jj. 



(28) 



If this is true, then a = —X/% + (y — fx)/2 when v > 0, and 
a = —1/6 — (3/2 when v < 0. Comparing with the locality 
condition in Eq. (f2"Tj) , we see that it is always satisfied. 
Based on numerical evidence we therefore conclude that 
the spatially extended system is able to adapt itself to 
variations in the exponents /i and v so that the locality 
criterion is always satisfied. 



D. Lack of dependence on spatial transport 
mechanism 

An important prediction of CFR is the lack of depen- 
dence on the diffusion constant. In Fig. [S] we show two 
sets of data for the same kernel but different diffusion 
constants. In one the diffusion constant is independent 
of mass. In the other it goes as D(m) tx m" 1 ! 2 such 
that k = —1/2. As can be seen, ir^m) scales exactly 
the same. However, as shown in the inset, the one point 
distribution function (N(m)) scales differently with k. 
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VI. HIGHER ORDER AGGREGATION 
PROCESSES 

Higher order aggregation processes may be considered 
where coalescence can only occur when n — 1 > 2 parti- 
cles meet at a single site. Although such processes have 
fewer physical applications than the binary case (n = 3), 
they have been suggested as an appropriate model of cer- 
tain polymeric reactions [20j and have received some at- 



J 



tention in the literature [ljj, HH [HJ . From a theoretical 
perspective, such systems provide an illustrative example 
of the breakdown of the CFR scaling due to a violation 
of the locality criterion. For these reasons, we consider 
the extension of the CFR argument to such systems. 

Again, we will restrict ourselves to the case where the 
reaction rate A(roi, . . . , m n _i) is a homogeneous func- 
tion of its arguments of degree (3. The Hopf equation 
corresponding to Eq.([5]) is: 



dt 



)/> CO ^ 
N(m) = / T\dm i X(m 1 ,...,m n - 1 )T[N(m i )d 
J o i=l i=l 



-(n-1) / IT dm,A(TOi, . . . ,m n ^ 2 . m )N(m) TT A(m l 
70 i=i <=i 



) + — (S(m-mo), (29) 
m 



r 



The flux-carrying correlation function is the (n— l)-point 
correlation function denoted by 

C(mi, . . . ,m„_i) = (N(mi) . . . iV(m n _i)}, mi ^ m^-. 

(30) 

By analogy with Eq.®, we introduce a quantity 
/(mi, . . . ,m„_i;m„): 



/(mi, . . . ,m„_i;m„) = 
A(mi, . . . , m„_i)C(mi, . . . , m n _i)5 



^ mi - m„(31 



On taking average in Eq. (|29p . the diffusion term drops 
out. Then, for m > mo we can write Eq. (|29[) as 



d(N(m)) 



dt 

n-l 



IT 



/(mi, . . . ,m n _i;m) 



/(mi, . . .,mj-x,m,mj+i, . . .m n -i;mj) 



The Zakharov transformations are:. 

mrrii 



rri; 



in; 



i + 3, 



(32) 

(33) 
(34) 



one for each of the n — 1 negative integrals. They have 
Jacobians (m/mj) n . Looking for homogeneous solutions, 

C(Ami, Am 2 , . . . , Am n _i) = A' l C(mi, m 2 , . . . , m n -i) 

. (35) 

and using the homogeneity exponent of A, we obtain 



0=/ T\dm i I(rni,... ) m n -i;m) 

J ° i=l 



n-l 



(36) 



where y = —h—0—n+l. We obtain a stationary solution 
when h — —j3 — n. The uniqueness argument of Sec. IIIII 
is easily extended to the case of n-ary interactions. 

It is cumbersome to discuss in full generality locality 
for higher values of n. Instead, we do a mean field analy- 
sis for n — 4 (three particles coalesce to form a new par- 
ticle) for the additive kernel A(mi,m,2,m.3) = mf + mf + 
TOg. In the mean field limit (N (mi) N (1712) N (1113)) ~ 
(minima)' 1 / 3 . When considering the collision integrals 
as a function of mi (say when mi — > 00), there is a free 
integral over m 2 . This integral being an integral over a 
pure power law, will either diverge at 00 or at 0. Hence 
the integrals are no longer finite and the locality con- 
dition will not be satisfied. Physically what happens is 
that three body collisions between three large particles 
are overwhelmed by three body collisions involving two 
large particles and one particle of very small mass. Thus 
the system behaves effectively as n = 3. One can get over 
this problem by introducing local kernels as discussed be- 
low. 

We now present some numerical results for n = 4. We 
consider additive kernel with {3 = 0, i.e, 



A(mi,m 2 ,m 3 ) = 1. 



(37) 



and measure the quantity 



7r 3 (m) 



dm 1 dm 2 {N(m)N(m 1 )N(m 2 )}. (38) 



,-2-/3 



which has a constant flux scaling of tt^ (m) 

For P — 0, the upper critical dimension is one. Hence, 
by the argument above, the locality condition should be 
violated and we should get scaling corresponding to n = 
3. In Fig. O we show the variation of ^(m) with m. 
The bottom curve correspond to the above kernel. CFR 
predicts that Tr^m) ~ mT 2 The bottom curve scales as 
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FIG. 6: The variation of n^m) with m is shown for the non lo- 
cal kernel Eq. (|37[) [bottom curve] and the local kernel Eq. (|39|) 
[top curve]. The simulations are for n = 4. The solid lines 
have slope —3 and —2. Inset: the variation of (N(m)) with 
m is shown for the local and nonlocal kernels. Their scaling 
with m is independent of the kernel for large m. 



m corresponding to scaling as predicted by /3 = and 
n = 3. 

To restore CFR, we consider a local kernel of the form 

A(TOi,TO 2 ,m 3 ) = g ( — ) g I — ] g ( — ] , (39) 
\m 2 J \m 3 J \mij 

where the dimensionless function is chosen to be 

1 



g(x) — exp [ iH 2 



(40) 



This local kernel has the effect that it suppresses interac- 
tions between masses that are not of the same magnitude. 
The results of 7r3 (m) for this local kernel is presented in 
the top curve of Fig. [SJ As can be seen, CFR is now 
obeyed. The inset of Fig. [6] shows that for both the local 
and non-local kernels (N(m)) has the same scaling, this 
is again as expected because both for n = 3 and n = 4, 
(N(m)) ~ rn~ 4 / 3 modulo log corrections for n = 4. 



VII. 



SUMMARY AND CONCLUSIONS 



To summarise, we have performed an extensive theo- 
retical and numerical study of the applicability and con- 
sequences of the CFR argument introduced in [14( in the 
context of cluster-cluster aggregation with a monomer 
source. We have used a heuristic scaling argument and an 
exact anaylsis of the appropriate Hopf equation to show 
that the scaling of the flux-carrying correlation function 
in the stationary state is fixed by the fact that the ele- 
mentary coalescence interactions conserve mass. In the 
case of cluster-cluster aggregation, the flux carrying cor- 
relation function is proportional to the probability of n— 1 



clusters coming together at the same point in space. It 
is thus not an esoteric object but is of direct physical 
significance. 

The CFR scaling exponent is identical to that given by 
mean-field theory. It is thus independent of the physical 
dimension and independent of the details of the spatial 
transport mechanism. This latter fact we have demon- 
strated clearly with some numerical simulations of aggre- 
gation with mass-dependent diffusion rates. The impor- 
tance and non-triviality of the result lies in the fact that 
the flux-carrying correlation function exhibits the mean 
field scaling even in the diffusion limited regime where 
mean field theory fails to give correct answers for other 
correlation functions, in particular for the density. This 
runs counter to the usual intuition in interacting particle 
systems where it is canonical that statistics are domi- 
nated by diffusive fluctuations in low dimensions where 
mean field theory breaks down. We do not consider our 
result to be at odds with this canon. It is indeed the case 
that most statistical quantities measured in the diffusion 
limited regime will be fluctuation dominated. What we 
have shown is that there is a particular special correlation 
function which does not feel these fluctuations at all. 

The usefulness of this result has already been demon- 
strated in our earlier work [TJ [l^] on constant kernel 
aggregation in low dimensions where it allowed us, taken 
together with a known exact result for the density, to 
prove multiscaling for the statistics of constant kernel ag- 
gregation in one dimension. Given the very direct phys- 
ical meaning of the flux-carrying correlation function in 
the aggregation context, it seems likely that other ap- 
plications will arise in concrete problems. At the very 
least, one can envisage using the result as a benchmark 
for numerical simulations of more complicated aggrega- 
tion problems, much as the 4/5-Law is used in validating 
numerical simulations of turbulence. 

As we have stated in our earlier paper, the CFR is not 
a theorem. It requires that a criterion which we refer 
to as "locality" should hold. To reiterate, by "locality", 
we do not mean that only clusters of equal masses are 
allowed to coalesce (although such a restriction would 
certainly ensure that our criterion is satisfied). Rather 
we mean a much weaker requirement that the mass in- 
tegrals describing the flux should not be dominated by 
their upper or lower limits. In general, locality is not 
testable a-priori. We have therefore devoted a consid- 
erable amount of effort in this article to studying the 
locality criterion in the context of cluster-cluster aggre- 
gation. From the theoretical perspective, we showed that 
if scaling exponents describing a system satisfy locality 
at the mean-field level (something which can be checked 
apriori) then there is a perturbative neighbourhood of 
models below the critical dimension for which locality 
holds. We then showed numerically that it is satisfied for 
a range of kernels in one dimension but breaks down for 
kernels for which one would expect long-range (in mass 
space) interactions to become dominant. We provided an 
instructive illustrative example of how the breakdown in 
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locality may violate CFR using a model kernel where the 
long range interactions may be tuned. 

It is rare that a generic nonequilibrium system will be 
solvable as the model discussed in the paper. It could 
be that the distinction between driving and dissipation 
scales get fuzzy (23|, or it could be that identifying the 
conserved quantity is a problem. In a recent paper [24| . 



we studied a model wherein the dissipation scale in not 
very well defined, and conjectured a CFR for such a 
model, even though it would not be expected apriori. The 
consequences of this conjecture was verified numerically. 
It would be of interest to clarify these observations the- 
oretically so that the results of the present article might 
be extended to an even wider class of models. 
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